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Abstract 

00 ' 

The adoption of two distinct boundary conditions for two fermions species on a finite 
lattice allows to deal with arbitrary relative momentum between the two particle 
CN . species, in spite of the momentum quantization rule due to a limited physical box 

size. We test the physical significance of this topological momentum by checking 
in the continuum limit the validity of the expected energy-momentum dispersion 
■ relations. 
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1 Introduction 



Among the restrictions of field theory formulations on a lattice, the finite 
volume momentum quantization represents a severe limitation in various phe- 
nomenological applications. For example, in a two body hadron decay where 
the energies of the decay products, related by 4-momentum conservation to 
the masses of the particles involved, cannot assume their physical values un- 
less these masses are consistent with the momentum quantization rule. In 
this letter we propose a solution to the problem based on the use of different 
boundary conditions for different fermion species 1 . 

We test the idea in the simplest case of a flavoured quark-antiquark cor- 
relation used to determine asymptotically the energy of the corresponding 
meson. In this case the fermion and the antifermion are the different fermion 
species and we show that suitable different boundary conditions can propagate 
a meson with a momentum that can assume continuous values. 

Section 2 introduces the boundary conditions, section 3 reports on the 
numerical results and section 4 summarizes the conclusions. 



R.P. thanks M. Liischer for drawing his attention on this point. 



2 Generalized boundary conditions 



In order to explain the method to have continuous physical momenta on a 
finite volume we first re-derive, for the sake of clarity, the momentum quanti- 
zation rule in the case of a particle with periodic boundary conditions (PBC). 
To this end we consider a fermionic field ip(x) on a 4-dimensional finite volume 
of topology T x L 3 with PBC in the spatial directions 

ip(x + Ci L) = ip(x) , i = 1,2,3 (1) 

This condition can be re-expressed by Fourier transforming both members of 
the previous equation 

J d 4 p e~ ip ( x+Si L) ^(p) = J d 4 p e~ ipx ^(p) , i = 1,2,3 (2) 

It follows directly from the previous relation that, in the case of periodic 
boundary conditions, one has 

e^ L = l =► P, = ^, .= 1,2,3 (3) 

where the n^s are integer numbers. The authors of [1] have first considered a 
generalized set of boundary conditions, that here we call ^-boundary condi- 
tions (0-BC), depending upon the choice of a topological 3-vector 9 

ip(x + 6i L) = e i6i ip(x) , i = 1,2,3 (4) 



The modification of the boundary conditions affects the zero of the momentum 
quantization rule. Indeed, by re-expressing equation (4) in Fourier space, as 
already done in the case of PBC in equation (2), one has 

«*.-*>* = 1 => P, = £ + ^, 4 = 1,2,3 (5) 



It comes out that the spatial momenta are still quantized as for PBC but 
shifted by an arbitrary continuous amount (Oi/L). The observation that this 
continuous shift in the allowed momenta it is physical and can be thus prof- 
itably used in phenomenological applications is the key point of the present 
work. The generalized ^-dependent boundary conditions of equation (4) can 
be implemented by making a unitary Abelian transformation on the fields 
satisfying #-BC 

idx 

ip(x) — >■ U (0,x)tp(x) — e~~ ip(x) (6) 



As a consequence of this transformation the resulting field satisfies periodic 
boundary conditions but obeys a modified Dirac equation 
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^(x) D e (x,y) ^(y) (7) 
x,y 

where the ^-dependent lattice Dirac operator D e (x,y) is obtained by starting 
from the preferred discretization of the Dirac operator and by modifying the 
definition of the covariant lattice derivatives, i.e. by passing from the standard 
forward and backward derivatives: 

V^(x) = - [U^(x)ijj(x + a /}) - ip(x)] 

CL 

Vpif)(x) = - tp(x) — U^(x — a p)i/)(x — a ft) (8) 

(X 

to the ^-dependent ones 

V n{e)il){x) = - [A M U^x)i>(x + ap,)- i){x)\ 

(X 

V M (0)ty(a;) = - [*l>(x) - K 1 U ^( x ~ a Mi* - a A)] (9) 

CI 

where we have introduced 

K = e— #o = (10) 

The authors of ref. [2] have considered for the first time #-BC in pertur- 
bative phenomenological applications. They used the shift in the momentum 
quantization rule, that they called a "finite size momentum", in order to build 
an external source to probe the tensor structure of the Wilson operators. 
A similar analysis was then repeated non-perturbatively by the same group 
in ref. [3]. The use of 0-BC has been considered in different contexts also 
in [4,5,6,7,8]. 

— * 

In this work we point out that the term 9/L acts as a true physical mo- 
mentum. 

As a test, we calculate the energy of a meson made up by two differ- 
ent quarks with different 8-BC for the two flavours. We work in the 0(a)- 
improved Wilson-Dirac lattice formulation of the QCD within the Schrodinger 
Functional formalism [9,10] but, we want to stress that the use of 8-BC in the 
spatial directions is completely decoupled from the choice of time boundary 
conditions and can be profitably used outside the Schrodinger Functional for- 
malism, for example in the case of standard periodic time boundary conditions. 
Let us consider the following correlators 

fm xo) = - y £ < CiW-YsCiW tjfrhMx) ) (11) 

y,z,x 
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where % and j are flavour indices, all the fields satisfy periodic boundary con- 
ditions and the two flavours obey different ^-modified Dirac equations, as 
explained in equations (7), (8) and (9). In practice it is adequate to choose 
the flavour i with = 0, i.e. with ordinary PBC, and the flavour j with 9^0. 
After the Wick contractions the pseudoscalar correlator of equation (11) reads 

fi(0;x o ) = y £ Tr ( 75 S 3 (6;z,x) 7s S t (0;x,y) ) (12) 

y,z,x 



where S(9] x, y) and 5(0; x, y) are the inverse of the ^-modified and of the 
standard Wilson-Dirac operators respectively. Note that the projection on 
the momentum 9/L of one of the quark legs in equation (12) it is not realized 
by summing on the lattice points with an exponential factor but it is encoded 
in the ^-dependence of the modified Wilson-Dirac operator and, consequently, 
of its inverse S(9; x, y). 

This correlation is expected to decay exponentially at large times as 

fA0;X O ) ^ f.. e -axoE ij{ e,a) (13) 



where, a part from corrections proportional to the square of the lattice spacing, 
Eij is the physical energy of the mesonic state 



Eij{6,a) = 



M* + [ 9 T \ +0(a 2 ) 



(14) 



here is the mass of the pseudoscalar meson made of a % and a j quark 
anti-quark pair. In the next section we will show the calculation of the meson 
energies for different flavours and for different choices of 9. We will show 
that after the continuum extrapolations we will find the expected relativistic 
dispersion relations 

4 = Ml + (i\ (15) 



3 Numerical tests 



All the results of this section are obtained in the quenched approximation 
of the QCD. We have done simulations on a physical volume of topology Tx L 3 
with T = 2L and linear extension L = 3.2 r , where r is a phenomenological 
distance parameter related to the static quark anti-quark potential [11]. In 
order to extrapolate our numerical results to the continuum limit we have 
simulated the same physical volume using three different discretizations with 
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0.362(5) 
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0.134208 


0.655(9) 


6.211 


24 


0.134540 


0.521(7) 






0.134954 


0.354(6) 






0.135209 


0.251(5) 






0.134517 


0.676(15) 


6.420 


32 


0.134764 


0.540(12) 






0.135082 


0.365(10) 






0.135269 


0.262(9) 



[0x,0 y ,0z] = [0.0, 0.0, 0.0] [1.0, 1.0, 1.0] [2.0, 2.0, 2.0] [3.0, 3.0, 3.0] 



Table 1 

Parameters of the simulations. The values of the bare couplings has been chosen in 
order to fix the extension of the physical volume L = 3.2 r$. For each value of the 
k parameter we have simulated all the values of 6. 



number of points (32 x 16 3 ), (48 x 24 3 ) and (64 x 32 3 ) respectively. We have 
fixed the three values of the bare couplings corresponding to the different 
discretizations using the r scale with the numerical results given in [12]. All 
the parameters of the simulations are given in table 1. The values of the RGI 
quark masses reported in table 1 have been calculated starting from the PCAC 
relation 

^pcac _ 

mu ~ 2f$(0;x ) {lb) 

where <9 M , M are the usual forward and backward lattice derivatives respec- 
tively while <9 M is defined as (<9 M + The time correlator fp(0;x ) has 
already been defined in equation (11) while /^(0;xo) is defined in the follow- 
ing relation 

a 6 - 

fi(9;x ) = -— ( Ci(y)l5(j(z) ipj{x)-fo-f 5 ipi(x) ) (17) 

y,z,x 

The improvement coefficient ca has been computed non-perturbatively in [13]. 
The RGI quark masses are connected to the PCAC masses of equation (16) 
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from the following relation 



,.RGJ 



m;~- = Z M (g ) [1 + (b A - b P ) ami] m^ AC (g Q ) 



(18) 



where the renormalization factor Zm(<?o) has been computed non-perturbatively 
in [14]. Also the difference of the improvement coefficients b A and bp is known 
non-perturbatively from [15,16]. In (18) the masses mi are the bare ones de- 
fined as 



arrij 



1 1 

ki kr 



(19) 



For each value of the simulated quark masses reported in table 1 we have 

— * 

inverted the Wilson-Dirac operator for three non-zero values of 9. Setting the 
lattice scale by using the physical value r = 0.5 fm, the expected values of the 
physical momenta associated with the choices of given in table 1 are simply 
calculated according to the following relation 



~ 0.125 101 GeV 



0.000 
0.217 
0.433 
0.650 



GeV L ~ 1.6 fm (20) 



These values have to be compared with the value of the lowest physical momen- 
tum allowed on this finite volume in the case of periodic boundary conditions, 
i.e. \p\ ~ 0.785 GeV. 

At fixed cut-off, for each combination of flavour indices and for each value 
of 9 reported in table 1 we have extracted the effective energy from the corre- 
lations of eq. (11), fp(9;x ), as follows 



a El J ff (9,a;x ) 



1 

2 



log 



7]^o-i) ' 

J l A9;x + l), 



(21) 



In fig. 1 we show this quantity for the simulation performed at /3 — 6.211 
corresponding to r mf G/ = 0.655 and r mf G/ = 0.354, for each simulated 

— * — * 

value of 9. As can be seen the correlations with higher values of \9\ are al- 
ways greater than the corresponding ones with lower values of the physical 
momentum 



|0l| > |0 2 



E l e J ff (9 u a;x ) > El J ff (9 2 ,a;x ) 



(22) 



a feature that will be confirmed in the continuum limit. 
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Fig. 1. Effective energies E l Jff(6,a;xo), as denned in eq. (21) at fixed cut-off. The 
results correspond to the simulation done at /3 = 6.211 with tq mf = 0.655 and 
tq rri2 = 0.354. Similar figures could have been shown for other combinations of 
the simulated quark masses and for the other values of the bare coupling. 



In the continuum extrapolations we have fixed the physical values of the 
quark masses slightly interpolating the simulated sets of numerical results. 
Being interested in the ground state contribution to the correlation of eq. (11), 
we have averaged the effective energies over a ground state plateau of physical 
length depending upon the quark flavours. We call (9, a) the result of the 
average and in fig. 2 we show a typical continuum extrapolation of this quan- 
tity. Similar figures could have been shown for the other values of simulated 
quark masses. 

The continuum results verify very well the dispersion relations of equa- 
tion (15) as can be clearly seen from fig. 3 in which the square of E l i(9) for 
various combinations of the flavour indices is plotted versus the square of the 
physical momenta \9\/L. The plotted lines have not been fitted but have been 
obtained by using as intercepts the simulated meson masses and by fixing their 
angular coefficients to one. 



7 



3 



_ 




c 


« 




* 
















# 


... • 






















----- — s 


S 













9 ' " 






|9| = 0.00 kj: _ 








jej = 1.73 »: 








|9| = 3.46 »i 








1 6 1 = 5.20 m 



0.01 0.02 0.03 0.04 0.05 

Fig. 2. Continuum extrapolations of the plateau averaged effective energies E %3 {Q, a). 
The results correspond to the quark masses tq = 0.650 and r^m^ 1 = 0.350. 

Similar figures could have been shown for other combinations of the simulated quark 
masses. 
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Fig. 3. Continuum dispersion relations. The data correspond to different combina- 
tions of the simulated quark masses and reproduce very well the expected theoretical 
behavior, i.e. straight lines having as intercepts the meson masses and as angular 
coefficients one (see eq. 15). 
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4 Conclusions 



We have argued that the limitation represented by the finite volume mo- 
mentum quantization rule can be overcame by using different boundary con- 
ditions for different fermion species. 

We have supported this observation by calculating the relativistic disper- 
sion relations satisfied by a set of pseudoscalr mesons in the case of quenched 
lattice QCD. We have shown that the physical momentum carried by these 
particles can be varied continuously by enforcing different ^-boundary condi- 
tions (see eq. (4)) for the two quarks inside the mesons. 

The method proposed can be applied to study all the quantities of phe- 
nomenological interest that would benefit from the introduction of continuous 
physical momenta like, for example, weak matrix elements. The suggestion 
can be applied in quenched QCD also in the case of flavourless mesons while 
can be extended to full QCD in the flavoured case only. 
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